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The interiors of neutron stars contain matter at very high densities, in a state 
that differs greatly from those found in the early universe or achieved at ter- 
restrial experiments. 1 Matter in these conditions can only be probed through 
astrophysical observations that measure the mass and radius of neutron stars 
with sufficient precision. 2 Here I report for the first time a unique determina- 
tion of the mass and radius of the neutron star EXO 0748—676, which appears 
to rule out all the soft equations of state of neutron star matter. If this object 
is typical, then condensates and unconfined quarks do not exist in the centers 
of neutron stars. 

The neutron star source EXO 0748—676 has repeatedly shown thermonuclear X-ray 
bursts. It has also shown, first 3 with EXO SAT and very recently 4 with RXTE, a special 
type of thermonuclear X-ray burst that is strong enough to lift up the outer layers of the 
star. During these so-called photospheric radius expansion 5 bursts, the radiation flux that 
emerges from the stellar surface is limited by the Eddington flux. It was, however, the 
detection 6 of redshifted O and Fe lines by XMM-Newton from the surface of the neutron 
star that rendered this low-mass X-ray binary unique in its class. For the first time, 
multiple phenomena have been observed from a single source that can be concurrently 
used to determine uniquely the equation of state of its dense interior. 

The dependence of the three observables, i.e., the Eddington limit, the redshift, and 
the ratio F cool /T c 4 , which is a quantity closely related to the total emitting area, on the 
stellar properties is shown in Table 1. Even though each of these measurements involve 
a different combination of the neutron star mass, radius, and distance, when all three 
observables are known, the expressions can be combined and inverted to yield a unique 
solution for the three stellar parameters as shown in Table 2. This is also shown visually 
in Figure 1, where the different constraints on the neutron star mass and radius imposed 
by each observable intersect at a unique set of values. One additional observable that, 
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in principle, can tighten these constraints is the amount of rotational broadening in the 
redshifted lines. This leads to a direct measurement of the radius of the star, 7 with an 
uncertainty arising from the angle i that the observer makes with the rotational axis of 
the star (Fig. 1). 

The only model parameters that are required in this solution are the color correction 
factor and the electron scattering opacity n es . The color correction factor f^, for an 
observer at infinity, relates the color temperature to the effective temperature of the star 
by Zoo = T c /T c fj and can be obtained by theoretical modeling of bursting neutron-star 
spectra. We use here a fitting function we devised that describes to an accuracy of 4% 
the results of recent model atmosphere calculations 8 for solar-abundance composition 



Note that foo in Tables 1 and 2 corresponds to the color correction factor at the cooling 
tails of the X-ray bursts when the emerging flux is substantially sub-Eddington. At this 
limit, models of bursting atmospheres 8 show that asymptotes to a value of ~ 1.34, 
practically independent of composition, temperature, and stellar gravity, significantly 
reducing the model dependence of our results. For the electron scattering opacity we 
use k cs = 0.2(1 + X) cm 2 g _1 , where < X < 1 is the hydrogen mass fraction. This 
is appropriate for the fully ionized neutron star atmosphere during the radius-expansion 
episode. In this calculation, we allow for the entire range of values for X. 

The main systematic uncertainty in our results arises from the assumption that the 
thermonuclear flash engulfs the entire star during the radius expansion and cooling phases 
of the bursts. There are a number of theoretical and observational arguments that support 
this assumption. First, numerical models 9 of the spreading of thermonuclear bursts on 
the surface of rotating neutron stars show deflagration times that are << 1 s, which is 
shorter by orders of magnitude than the duration of the bursts. Second, the constraints 10 
on the magnetic field strength of EXO 0748—676, imposed by the lack of Zeeman splitting 
of the atomic lines, show that the magnetic field is dynamically unimportant and thus 
cannot inhibit the spreading of the thermonuclear flash over the entire surface. 

On the observational side, further evidence for uniform emission from the entire surface 
as well as the constancy of the Eddington flux during the bursts is obtained from the 
study 11 of the peak luminosity of a large number (~ 70) of bursts from 4U1728— 34, which 
shows that the peak flux is constant to within 1% in bursts separated by months. Indeed, a 
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larger study 12 of peak fluxes of all thermonuclear burst sources in the RXTE catalogue also 
show a quantitatively similar result. Second, in the cooling tails of thermonuclear bursts, 
the observed ratio F coo i/T 4 asymptotes to a constant and reproducible value between 
bursts, strongly indicating that the entire surface of the neutron star emits uniformly in 
the cooling tails. 13,14 

Even though flux oscillations of amplitude <10% have been seen in some thermonu- 
clear bursts and are thought to be caused by modes excited on the neutron star surface 
during the bursts, 15 their presence does not increase the uncertainties reported here. First, 
bursts oscillations have never been observed during the radius expansion phases of the 
bursts but only in the rise phase and cooling tails, thus not affecting the determination 
of the Eddington limit from observations. Second, the low amplitudes of the oscillations 
during the cooling tails introduce uncertainties that are smaller than both the system- 
atic and statistical uncertainties (10%) allowed for in our calculations. This is especially 
true for the burst oscillations in EXO 0748— 676, which are very weak, with an average 
amplitude of 3%. 16 

A final consideration about the Eddington limit is related to the point during a radius 
expansion burst at which this quantity is measured. The relevant measurement here is 
the so-called touchdown flux, which is the Eddington limit at the point in the burst when 
the photosphere returns to the actual radius of the star. 14 The flux at the touchdown 
point cannot be smaller than the Eddington flux because, if the radiation support of the 
atmosphere were suddenly removed, the photospheric radius would return to the actual 
size of the star within a free-fall timescale, which is ~ 1 ms for a neutron star. Instead, an 
adiabatic return to the stellar surface, at timescales of ~ few seconds, is observed in every 
radius expansion burst, indicating that the emerging flux traces the Eddington limit all 
the way to the touchdown point. 

The only unknown property of the binary system that affects the mass and radius 
measurements is the hydrogen fraction X of the accreting material. Further observations 
of the elemental abundances of the accretion flow at long wavelengths can greatly reduce 
this uncertainty. However, even when we consider the most extreme value of X = 0.7, we 
obtain lower limits on the mass and radius: 

M > 2.10 ± O.28M , R > 13.8 ± 1.8 km, (2) 

which is shown in Figure 2. For smaller values of the hydrogen mass fraction, i.e., 
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for a Helium-rich companion which can be expected in a small-orbit binary such as 
EXO 0748—676, the values of M and R are even higher. It is clear from this figure 
that the unknown value of X is by far the dominant systematic uncertainty in our result. 
Also note that even though separate uncertainties on mass and radius are quoted, these 
are not independent and do not form an error ellipse because the mass-to-radius ratio is 
fixed by the redshift measurement with a negligible error. 

Because we treat distance as an independent variable, we can use the same three 
measurements to determine a lower limit on the distance to the source, D > (9.2±1.0) kpc. 
This is a significant advantage of our method: The equations in Table 2 show that if the 
redshift is measured, the mass and radius measurements are independent of the distance 
and vice versa. Alternatively, if there is a direct measurement of the distance, e.g., for a 
source in a globular cluster, then only two spectroscopic measurements are sufficient to 
determine the mass and radius of the neutron star. 

The mass-radius measurement presented here does not suffer from uncertainties related 
to magnetic and centrifugal support that can complicate the calculation of neutron star 
structure given an equation of state because of the slow rotation and the low magnetic 
field of EXO 0748— 676. The large size and the high mass of this neutron star impose 
stringent constraints on the equation of state of matter at supernuclear densities (see 
Figure 2). Even the lowest allowed value for the mass and radius can be obtained by 
only the stiffest equations of state. In particular, equations of state with condensate 
interiors predict large radius-, low- mass neutron stars (Fig. 2) that are mostly excluded 
by our measurements. On the other extreme of the mass-radius diagram, self-bound 
bare strange-matter stars have mostly small masses and radii, also inconsistent with 
the mass and radius of EXO 0748— 676. The more conventional equations of state with 
neutron-proton compositions are more likely to explain the set of observational properties 
of EXO 0748—676 presented here. This argues that hadrons and not deconfined quarks 
represent the ground state of matter. 
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Table 1. The three main quantities observed from EXO 0748— 676 and their theoretical depen- 
dence on the neutron star properties. The Eddington limit -FEdd> defined as the radiation flux 
at which the outward radiation force balances the inward gravitational force, is the limiting flux 
emerging from thermonuclear X-ray bursts with photospheric radius expansion. The measure- 
ments of the touchdown flux reported here were obtained by averaging the values determined 
recently with RXTE A and earlier with EXOSAT 3 observations, which are consistent with each 
other. The redshift z of O and Fe absorption lines in the X-ray burst spectra of EXO 0748—676 
has been measured for the first time with XMM-Newton. 6 The ratio F coo \/aT A , where -F coo i 
and T c are the thermal flux and the color temperature inferred from the X-ray burst spectra, 
respectively, asymptotes to a constant value during the cooling tails of the bursts. This ratio is 
directly related to the emitting surface area of the neutron star and was inferred from the RXTE 
data 4 by taking the averages of F coo \/aT* measured between 7—15 s after the burst, ensuring 
that a constant ratio has been reached following the photospheric radius expansion but that the 
observed flux is still high enough for an accurate determination. The physical constants G and c 
are the gravitational constant and the speed of light, respectively, K es = 0.2(1 + X) cm 2 /g is the 
electron scattering opacity, X is the hydrogen mass fraction of the accreted material, is the 
color correction factor, and D is the distance to the source. All the above expressions include 
minimal general relativistic corrections that are accurate in the slow rotation limit, which is 
appropriate for EXO 0748— 676 given its 47 Hz spin frequency. 16 



Table 2. The neutron star properties that are obtained from the observations summarized 
in Table 1. The neutron star mass, radius, and distance are uniquely determined from the 
Eddington luminosity, the redshift, and the ratio F coo \/aT A as shown in the middle column, 
given a model for the neutron star atmosphere that determines the color correction factor and 
a measurement of the hydrogen mass fraction. The last column shows the minimum values for 
the mass, radius, and distance to the neutron star for any possible value of the color correction 
factor and the hydrogen mass fraction. The counter-intuitive dependence of the stellar properties 
on the observables arises from the combination of the expressions shown in Table 1. These 
expressions also show why a unique spectroscopic determination of the mass and radius of a 
neutron star, which is necessary to constrain its equation of state, can only be achieved when 
all three phenomena are observed from a single source, if the distance to the source is unknown. 
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Figure 1. Four complementary methods to determine the mass and radius of a 
neutron star. The contours on the mass-radius plane of neutron stars imposed by the mea- 
surement of (i) the Eddington flux during photospheric radius expansion bursts (blue), (ii) the 
ratio F coo \/T^ of the surface emission that asymptotes to a constant during the cooling tail of 
a thermonuclear burst (green), (Hi) the redshift of atomic absorption lines observed during the 
burst (red), and (iv) the broadening of such lines due to the rotation of the star (magenta) for 
a hypothetical star with M = 1.8M@ and R = 10 km. The second quantity, obtained from 
the thermal flux F coo i and the color temperature T c of the burst spectrum, is closely related to 
the total emitting area from the stellar surface when the nuclear burst has engulfed the entire 
star. The widths of the contours correspond to a hypothetical 10% uncertainty in each measure- 
ment. The uncertainty in the redshift measurement can be negligible if a grating spectrometer 
is used, as in the case of EXO 0748—676. The black shaded area is the intersection of the four 
contours and corresponds to the uncertainty in the measurement of the true mass and radius of 
the neutron star. 
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Figure 2. The constraints on the neutron star equations of state imposed by obser- 
vations of EXO 0748—676. The predicted mass-radius relations for a number of representative 
equations of state of neutron stars without condensates (blue), with condensates (magenta), and 
for strange stars (red). 17 ' 18 The curve labels and the corresponding references can be found in 
refs. [15] and [16]. The two 1 — a error bars corresponds to the measurement of the mass and 
radius of EXO 0748—676 from three independent observations (equation 2) for different values 
of the hydrogen mass fraction X of the accreted material. The measurement for X = 0.7 cor- 
responds to the minimum allowed values for the mass and radius of the neutron star. Only the 
stiffest equations of state are consistent with this measurement. 



